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Abstract 

Inspired from Sutherland's work [Phys. Rev. Lett. 74, 816 (1995)] on 
detecting bound spin waves, we propose that bound electron states can be 
detected from the dependence of interacting electron systems to the Aharonov- 
Bohm flux in the 'extended zone' scheme, where the electron pairing halves 
the original period N a flux quanta in a system of linear size N a . Along with 
the Bethe-ansatz analysis, a numerical implementation for keeping track of the 
adiabatic flow of energy levels is applied to the attractive/repulsive Hubbard 
models and the t — J ladder. 

PACS numbers: 74.25.-q, 03.65.Ge, 71.27.+a 
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The response to an adiabatic change against external parameters is an interesting way to 
probe the nature of interacting electron systems. Some thirty years ago, Byers and Yang |I[ 
proposed a notable example of detecting the Cooper pairing: while a normal state responds 
to an Aharonov-Bohm (AB) flux periodically with the period of flux quantum, $o = hc/e, 
a BCS state will have a halved period, $o/2. The anomalous flux quantization has actually 
been applied to various strongly correlated electron systems that are intended to describe 
high-T c cuprates 

This is in fact one out of several ways to detect superconductivity in purely electronic 
systems. Since an effective electron-electron attraction per se does not guarantee the Cooper 
pairing, a usual way is to search for a long-tailed pairing correlation function, but the finite- 
size effect must be carefully analyzed. The quantum Monte Carlo method (QMC) 0-||, 
the density-matrix renormalization group 0, evaluation of the superfluid density (helicity 
modulus) flO|fl are toward this line of approach. Thus the detection of the Cooper pairing 
(or bound electrons in more general terms) is a demanding problem in a correlated electron 
gas, except for one-dimensional(lD) systems, where exact analytic treatment is feasible with 
the Bethe-ansatz solutions coupled with the conformal field theory. 

This is where the Byers- Yang flux quantization comes in. The test, however, remains 
some way from a clear-cut criterion, since the flux dependence of the ground-state energy 
may give the half- flux periodicity even for the repulsive Hubbard model 0. An origin of 
the obscured period is the spin degrees of freedom ||11|| , which is most prominent in ID. 



In the present paper, we propose a new way to detect bound electron states from a more 
global look at the response to the flux. The idea has been inspired by a recent analysis of 
the bound spin waves by Sutherland [|1^]: bound complexes of spin waves in ID Heisenberg 
magnets (or equivalently a gas of charged bosons) may be detected from their response to 
a boosted total momentum. The momentum boost is achieved by twisting the boundary 
condition, \l/(. . . ,Xj + N a , . . .) = e* 2?r ^\l/(. . . , Xj, . . .), for an A a -site lattice, which uniformly 
shifts the set of k points. He has shown that, if all of the N particles (i.e., N flipped spins 
in a magnet) form one bound state, the energy returns to its initial value by a twist of 



= N a /N, which is 1/N times the twist 4> = N a required to shift the set of k points for the 
free particles back to the original position. Intuitively, this discerns whether the momentum 
boost acts on individual particles or on a 'center-of-mass' of bound particles. A key in 
Sutherland's idea is to keep track of the 0-dependence of the state not over the one period 
(0 < (j)< 1) but over the 'extended zone' (0 < < N a ). 

Sutherland does not argue what happens in electron systems or in the situation where 
more than one bound complexes exist, but we conjecture here that the extended AB method 
should, conceptually, hold in general. We would then be able to treat e.g. the Cooper pairing 
problem. Natural questions are: (i) Can we really extend Sutherland's spin-wave analysis 
to electron systems? (ii) Can we apply the method to two or higher dimensions? 

In the present Letter we first give a straightforward extension of Sutherland's spin anal- 
ysis to electron systems with the Bethe-ansatz analysis of the ID Hubbard model as a 
prototype. There we introduce an AB flux $ that couples to the charge degrees of free- 
dom to twist the boundary condition. We then look at the energy levels against $ over 
< $ < N a $> Q (which we call the 'extended AB' spectral flow) to discriminate the bound 
states, as opposed to the conventional wisdom that one period of $ suffices. We next 
propose a method to numerically implement the extended AB test for arbitrary systems 
including 2D systems to go beyond the Bethe-ansatz analysis. 

We start with confirming for the ID Hubbard model that the electrons have a reduced 
period of N a § /N as well for iV-bound states with the Bethe-ansatz analysis fL3 |. Consider 



an iV a -site ring (with N a even for simplicity) containing N electrons, and thread a magnetic 
flux $. A change of $ by $o shifts the k points exactly by their spacing Ak = 2ir/N a 
for noninteracting electrons, so that the set of k points accomplishes a full travel across 
the Brillouin zone when $/$o reaches N a . When interacting, the electrons are subject to 
Bethe-ansatz equations for which there are two types of solutions, i.e., real and complex 
roots ||15|| . A real charge rapidity kj represents the quasi- momentum of a beam of charges, 
while complex fc^'s, which are sometimes called string solutions because they appear in a 
linear group with a common real part, represent bound states of electrons. 



For a set of real kj's, we can take a logarithm of Bethe-ansatz equations to have 

M 4 

kjN a = 2nlj + 27T$/$ - ^ 2 tan" 1 — (sin kj - A a ) , (1) 

a=l ^ 

where M is the number of down spins, t the transfer, U the Hubbard interaction, and Ij 
(j = 1, . . . , N) is an integer (half odd integer) for an even (odd) M. We assume that spin 
rapidities A Q are real as is the case with the ground state of the repulsive model. Equation 
(HI) has a periodicity of N a in Ij. Since an increase of $ by $o causes a uniform shift of Ij by 
unity, a real solution indeed has a periodicity of (except at the half-filling, see below). 

Now, the string solutions for charge rapidities, {k n i (/ = l,...,2n)}, are specified by 
a single real parameter, A' n ||15|| . When all the TV particles form a single bound state (an 
./V-bound state in Sutherland's words), the equation reduces to 

exp (iN a Y^ Ki) = 
i 

exp {-iiVjsin- 1 ^ + inU/A) + sin" 1 ^ - inU/A))} = exp(i2nN$/<5> ) . (2) 

We can then see that is determined by the total momentum J2i k n i, so that a change of 
$ by &o/N is enough to shift A' n to the next position among the A^ a solutions. Thus we 
end up with a period of N a § Q /N for a single A^-bound state. 

Solutions having more than one set of strings have been known to exist for e.g., the 
attractive Hubbard model having a set of two-strings (electron pairs) [ |TB| , [TT| | . However, 
their high-energy spectra has not been fully understood, so that the spectral flow has to be 
obtained numerically even in ID. 

So we move on to the numerical implementation of looking at the flow, which is readily 
applicable to two or higher spatial dimensions, and also to various models such as the 
extended Hubbard or t-J models. In determining the extended period, we have to keep 
track of the ground state for the range of the flux well beyond the flux quantum, where 
the energy level soars to become high-energy states. In addition the level has to be traced 
straight through level crossings, so that the numerical method must be good enough to 
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reproduce (a) level crossings and (b) high-energy states. Conventional methods such as the 
exact diagonalization or quantum Monte Carlo methods would then be inadequate. 

The algorithm we propose consists of successive estimations of the energy and wave 
function for $ + A$. The new state * is estimated by multiplying a connection, e iA (*) A<J> ) 
defined as 

^($ + A$) = e iA ^ A **($). (3) 

Here a matrix A is approximately diagonal (e lAiA * ~ 1 £j($ + A$)/ , I'i( ( l ) )) for a small enough 
A$, where denotes the i-th component. Namely the connection reduces to a set of pure 
phases (i.e., Berry's geometrical phase). Each estimation is followed by a correction in the 
following recursive manner. 

i) Start from the ground state = 0) and the next one ^(A^) which are obtained by 
the Lanczos method. 

ii) Predict the connection from the previous step as e Jj4lA * ~ \I/j(<l>)/\l/i($ — A$) to produce 
the next one, + A$) = e tA ^ A *^($). 

hi) Correct \I/ by operating some power, m, of the resolvent as 

[t-, )— s 777 1 W + A$) T ™ + A$) (4) 

\ [H{§ + A$) — E] 2 J V ; v; 

This power method is performed with the inverse-iteration and the conjugate-gradient op- 
timization. E in the resolvent can be calculated accurately from the wavefunction updated 
every time we increase m. 

iv) Check the convergence, \(H(Q + A$) — E)^(Q + A$)| < e, for a given accuracy e. 
Increase m until this inequality is fulfilled, then move on to ii). 

Applied over the extended zone, the method can count the periodicity (or the 'winding 
number') of the state. The tractable matrix size is the similar to that for the conventional 
Lanczos diagonalization. 

Around level crossings the calculation becomes subtle. We first observe that level cross- 
ings come in two classes: One class occurs for two levels that differ in some symmetry of 
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the system. The second occurs at a phase transition point, e.g., the normal-superconductor 
transition. The first class can be dealt with by only letting the convergence criterion in iv) 
with typically e = 10~ 8 for the systems considered here with the matrix size up to a few 
ten thousands. The second class requires scaling analysis by varying the distance from the 
critical point as illustrated for the t — J ladder below. 

Our first example is the ID Hubbard model. Figure 1 compares the spectral flow for a 
10-site ring with 6 electrons for a repulsive model with U/t = +1 and an attractive model 
with U/t = —1. The extended period, which is N a $o for the repulsive model, is indeed 
halved for the attractive case, providing a clear-cut detection of the two-electron pairing. 

For comparison, we have superposed the anomalous flux quantization test for < $ < $o 
in Fig.l. We can see that the latter test is indeed ambiguous. To be more precise, we can 
observe the following. If the Fermi sea has a closed shell (N = 2 (mod 4), as is the case with 
Fig.l), the true ground state is always spin singlet irrespective of the sign of U. In addition 
to the branch starting from the ground state, there is a second branch, which is degenerate 
with the first one at $ = $o/2 for U = but shifts upward for a repulsive interaction, or 
downward for an attractive interaction. Thus a dip appears in a continuous fashion as a 
negative U is turned on. However, the half periodicity can appear even for the repulsive 
interactions for open-shell Fermi seas ||. This is due to the existence of a spin-triplet state, 
which is degenerate with the ground state at $ = and stabilized over the singlet state for 
U > fT8| . If we further go into the strong correlation limit, even more anomalous 1/N e - 



periodicity can appear as shown by Kusmartsev ||11|| . Thus we have to worry about these 
finite-size effects due to other branches lying around to obscure the period in the anomalous 
flux quantization. In two or higher dimensions such situations may be improved fL0|,|5||) 



but we should re-emphasize that the winding- number counting here concentrates on the 
adiabatic evolution of a single state, where the 'global' period is determined independently 
of other branches. 

The abrupt change from N a § to N a § /2 in the extended AB period takes place exactly 
at U = 0, which is the critical point in ID. The change occurs despite the fact that three 
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spectral flows (for U < 0, U = 0, U > 0) are almost identical around the respective minimum 



except for some offsets. This corresponds to the known fact [HJ that the charge stiffness (or 
Drude weight) does not exhibit singular jumps even at the critical point when the system is 
finite. 

A closer inspection shows that the long AB period for U ^ is dominated by some 
level crossings that turn into a level repulsion or level anti crossings, where different sets of 
anticrossings are selected according as U > or U < 0. It is at first puzzling how such 
a qualitative change can possibly occur for an infinitesimal \U\, since charge rapidity k^s 
coalesce into doubly occupied k points in the Fermi sea no matter how the critical point 
(U = 0) is approached from repulsive or attractive sides (see Ref. fL9f). 

The puzzle is resolved by looking at the weak-coupling Bethe-ansatz behavior. As the 
k points shift with $, the crucial level anticrossing occurs when the upper-most doubly- 
occupied k point reaches k = tt/2 situated at e(k ) = 0, the center of the band dispersion. 
This special point accommodates highly degenerate states that are connected by two-particle 
scattering processes that comprise the normal ones, <4 T c fcoJ. — > c L+pt c L-p!' w ^h p the 
momentum transfer, along with the Umklapp process, <4oT c fco4. — * c l-2fc T c fc-2fc J.' When U 
is switched on, the degenerate perturbation dictates that some of these states must mix to 
give anticrossings. 

The normal process with p = 2irN/N a selectively produces the anticrossing for the re- 
pulsive case. In contrast, the anticrossing in the attractive case is caused by the Umklapp 
process that transfers a pair to other e = points, which is the only possible process indeed 
if the pairs are not dissociated by the adiabatic change in $. The fact that ko is the key 
position is illustrated in Fig.l: it takes $ c /$o = (k — kF)/Ak for the upper-most doubly oc- 
cupied k point (kp) to reach ko, where $ c /$o = 1-5 for N e = 6, N = 10, in exact agreement 
with Fig.l. 

Given this property, we can readily show that a series of dissociations of doubly occupied 
states for U > gives the extended period of N a $ , while a series of Umklapped doubly 
occupied states for U < halves it. Thus, an infinitesimal interaction is enough to change 



the global topology (the winding number) of the connection in this example in ID. 

A fuller understanding in terms of the spin-charge separation in ID emerges if we look 
more closely at the Bethe ansatz, where A a (a 'spin degrees of freedom'), not directly coupled 
with $, stays approximately constant, while the 'charge degrees of freedom' ki progressively 
changes with $ for a repulsive U. This is accomplished by A a sequentially parting company 
with one ki to meet another, which is exactly where the anticrossings occur. On the other 
hand, each A a is attached to a pair, kj,k*, for an attractive U, and the center-of-mass 
momentum, 2Ke(kj), of the pair remains within the first Brillouin zone due to the Umklapp 
process, or equivalently, k/s have to satisfy cos(kj) > for the two-bound solutions to exist. 
The details will be published elsewhere. 

In order to demonstrate that the present method is applicable to the case in which the 
Bethe-ansatz is inapplicable, and also to provide a step toward higher dimensions, we move 
on to the t-J ladder model. This model, originally conceived for some copper or vanadium 
oxides, is believed to exhibit superconductivity for small hole doping from the half-filling 
ffijfl . Here we perform the extended AB test for a 6 x 2-site system with 4 electrons, which 



corresponds to a hole concentration of n = 2/3. 

In the result, Fig. 2, we can clearly see that the periodicity is halved into N a /2 as J 
is increased to J > 2.04£, which indicates paring. Here we can illustrate a nice feature 
about the extended AB test: since the transition is associated with a level crossing (a cusp) 
turning into an anticrossing, we can numerically plot (inset of Fig.2) the size, A, of the level 
repulsion against the relevant parameter (J here) to identify the critical point (J c ~ 2.04 
here) at which A vanishes. Thus we can estimate the critical value J c for a finite system in 
a well-defined manner, which may then be cast into a finite-size scaling analysis for a more 
pertinent definition of the critical point. Thus the present method provides a possible way 
to determine the phase diagram of a given model. 

So far we have not discussed anything directly about the coherence of the pairs, so that 
we are talking about a necessary condition for superconductivity. For that matter, Byers- 
Yang theorem also gives a necessary condition for the superconductivity (or the Meissner 
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effect). It is an intriguing future problem to see if the coherence can possibly appear in the 
spectral flow. 

Another comment is that the present test gives information on metal-insulator transitions 
such as the Mott-Hubbard transition as well. For the half-filled Hubbard model (a Mott's 
insulator), the extended AB period reduces down to $o- The sudden change from the full 
period N a <&o to the minimum one is due to an appearance of the charge gap, across which 
the flow is inhibited to jump so that the system has to return to the original state as soon as 
the flux reaches $o> i n analogy with the gauge argument of the quantum Hall effect. Thus 
we can expect that the present method can detect the existence or otherwise of the Fermi 
surface |[21|| . A more interesting example is the phase separation, where the system should 
respond as a single bound state as will be reported elsewhere. The different winding numbers 
for paired, metallic and insulating phases may be analyzed in terms of the homotopy of the 
phase-space (or the fiber bundle) of correlated electron systems. 
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FIGURES 

FIG. 1. The result for the spectral flow for a 10-site Hubbard ring with 6 electrons. The flow 
lines represent a repulsive case (U = t, solid line), an attractive case (U = —t, dashed line), and 
the non-interacting case (dotted line). The relevant level crossing point is marked with <& c . For 
comparison, the low-lying flows over one period of the flux quantum is reproduced. 

FIG. 2. The spectral flow lines for a t-J ladder model with 4 electrons in 6 x 2 sites for 
J = 2.2, 2.6, . . . , 3.8 ^from top to bottom. The inset depicts the size, A, of the relevant level 
repulsion against J, where the line is a guide to the eye. 
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